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Abstract 

We suggest a system in which the amplitude of macroscopic flux tunneling 
can be modulated via the Aharonov-Casher effect. The system is an rf-SQUID 
with the Josephson junction replaced by a Bloch transistor - two junctions 
separated by a small superconducting island on which the charge can be in- 
duced by an external gate voltage. When the Josephson coupling energies of 
the junctions are equal and the induced charge is q = e, destructive inter- 
ference between tunneling paths brings the flux tunneling rate to zero. The 
device may also be useful as a qubit for quantum computation. 
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Geometric-phase effects are ubiquitous in physics (T|]. The best known of these is the 
Aharonov-Bohm effect H in which charges moving in a field-free region surrounding a mag- 
netic flux nevertheless pick up a phase that leads to interference. This effect long predates 
the general formulation of the geometric phase by Berry [[|. A related effect was proposed 
by Aharonov and Casher [|J , who showed that a magnetic moment that moves around a line 
charge acquires a phase proportional to the linear charge density. Extending this analysis to 
superconductors, Reznik and Aharonov [0] showed that magnetic vortices (or fluxons) mov- 
ing around an external charge q could acquire an Aharonov- Casher phase. Other theoretical 
work has concerned the appearance of this effect in Josephson-j unction arrays and devices 
P^TTf. Experimental evidence for the Aharonov-Casher effect in a Josephson system was 
first provided by Elion et al. fl2| . There the relevant charge around which the vortices move 
is that induced on a superconducting island by an external gate voltage. Manifestations of 
the Aharonov-Casher effect in Josephson-j unction systems studied up to now - 2e-periodic 
modulation of an array resistance with the induced charge and quantization of the voltage 
across the array - are, however, not specific to this effect. Qualitatively similar phenomena 
can be produced simply by the Coulomb-blockade effect due to the quantization of charge on 



the superconducting islands of the devices [13]. In this Letter, we suggest a system, where 



the Aharonov-Casher effect manifests itself in a way that can be directly distinguished from 
Coulomb-blockade oscillations. 

The proposed system (Fig. |l]) is a variant of a SQUID. Macroscopic quantum tunneling 
of the magnetic flux in SQUID systems has been the subject of extensive experimental 
investigation in the past several years [|T4,[T5|] . Recent experiments have shown that a SQUID 
can be put into a coherent superposition of macroscopically distinct flux [IS] and persistent- 
current [ 16 1 states. The system we consider is a SQUID with the Josephson junction replaced 
by a Bloch transistor (see, e.g., JT3]) - two small, closely-spaced Josephson junctions, in 
which the charge on the island between the junctions can be varied by an external gate 
voltage V g . Because of the gate- voltage dependence of the critical current of the Bloch 
transistor, the flux-tunneling rate in such a SQUID should exhibit 2e-periodic oscillations as 
a function of the gate-voltage-induced charge q = C g V g . As we show below, in addition to 
these ("Coulomb-blockade") oscillations, the flux tunneling is completely suppressed when 
the induced charge is e. This fact does not follow from the simple charge quantization on 
the middle electrode of the Bloch transistor, but can be explained naturally in terms of the 
Aharonov-Casher effect. 

In a handwaving sense, flux tunneling occurs via the passage of a fluxon through one of 
the junctions and interference can occur between the two paths that encircle the induced 
charge q on the middle electrode. When q = e, the relative phase of the two paths is 7r and 
the interference is destructive. If, in addition, the Josephson coupling energies Eji&iad Ej 2 
of the two junctions (which control the amplitude of flux tunneling through each individual 
junction) are equal, Eji = E J2 fllT] , the destructive interference is complete, and the total 
amplitude of flux tunneling is exactly zero. More generally, the Aharonov-Casher effect in 
this system makes it possible to control the phase of the flux-tunneling amplitude, an effect 



that can simplify design of Josephson-j unction qubits for quantum computing [ |15| , |16| ,|i8H21 

Quantitatively, we consider a SQUID of inductance L (Fig. [I]) that is assumed to be not 
too large, so that the SQUID has only two metastable flux states p|J. The SQUID has two 



junctions with Josephson phase differences 4>\ and 2 across them. A gate electrode couples 
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with capacitance C g to the island between the junctions and an external flux <& x is applied 
to the SQUID loop. The Hamiltonian for this system is 

TT Q 2 (2en-q) 2 ($ - § x ) 2 „ 

where n is the number of Cooper pairs charging the island, Ce is its total capacitance relative 
to all other electrodes. The flux $ in the SQUID is related to the phase differences across 
the junctions by 27r$/$ — 4>\ + 02, where $ — h/2e is the flux quantum. $ is conjugate 
to the charge Q on the capacitance C between the ends of the SQUID loop: [&,Q] = ih. 
(Because of the stray-capacitance contribution to C, it is not related directly to the island 
capacitance Cj>) Similarly, the phase 9 = (<pi — $2)/ 2 is conjugate to n: [9, n] = i. Thus, $ 
and 9 are the "coordinates" of the SQUID and Q and n are the "momenta". 

One can then view the first two terms in the Hamiltonian as representing the kinetic 
energy of the system and the remaining terms as the potential energy, shown in Fig. |2] for 
the case $ 2 = $o/2 and Eji = Ej 2 . The potential has minima in a checkerboard-like pattern 
with minima on one side of $ = $o/2 shifted by 7r in the 9 direction relative to minima on 
the other side. The minima are not separated by integer multiples of $0 i n the $ direction 
because of the inductive term in Eq. (ffl). On the other hand, the potential is strictly periodic 
in the 9 direction and, consequently, neighboring minima are always separated by A9 = ±7r. 
In fact, since the number n of Cooper pairs on the island is integer, one can impose periodic 
boundary conditions on 9, wrapping all points with 9 = n to 9 = —n. 

We will analyze tunneling between different flux states in two limits, one in which the 
charging energy for the island Eq = (2e) 2 /2C*s is much smaller than Ej and the other 
in which Eq Ej. First, we will consider the former case, where 9 is a good quantum 
number and one can employ a tight-binding picture in which the system is localized in one 
of the minima of Fig. |^. The tunnel splitting can be found using an instanton technique in 
which the imaginary-time action Sj is evaluated along the least-action paths in the inverted 
potential [p3[| . The zero-temperature tunnel splitting is then given by 



A = 5> ie -**/\ (2) 
j 

where the sum is over all least-action paths between two potential minima, and Uj and Sj 
are the attempt frequency and imaginary-time action, respectively, for the jth path. The 
action is calculated by evaluating 

Sj = f C E (r)dT (3) 



along the jth path, where £e{ t ) = —C(t — ► — it) is the imaginary-time (Euclidean) La- 
grangian obtained from the real-time Lagrangian C, which in turn is obtained from the 
Hamiltonian by the usual transformation: £ = $Q + (|^#)(2en) — H. For Ej\ = Ej 2 = Ej 



lO (11)1 mil M' = (./ / ( ;uiu 1 " 

time Lagrangian is found to be 



and using Hamilton's equations to obtain $ = Q/C and -^9 = {2en — q)/C^ the imaginary- 



C E (r) = ^ 2 + ^[^) + ($ J x)2 - 2Ej cos(7r$/$ ) cos 9 - iq f , (4) 
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where the dot represents differentiation with respect to r. The last term in Eq. (P, being a 
total time derivative, can have no effect on the classical dynamics of the system. Yet, it has 
profound effects on tunneling, giving rise to the Aharonov-Casher phase and the resulting 
interference effect. 

To see this, let us consider the two tunneling paths illustrated schematically by the 
arrows in Fig. |2|. The endpoints of the two paths are equivalent, since they merely differ in 
9 by 27T. Every term in Eq. (0) except the last is symmetric under the reflection 9 — > —9. 
Hence, the actions of the two paths differ only because of the last term. The imaginary-time 
action for each path can then be separated into two parts: Sj = Si + S J geo , where Sj is the 
action obtained by using all but the last term in Eq. and is the same for both paths. 
The geometric-phase action for each path is given by 

S]i = ~n (fj) / Odr = T iq = T™(q/2e)h. (5) 

1,2 

From Eq. (0), the tunnel splitting is given by 

A = 2A cos(g7r/2e), (6) 

where A = cuoe~ Sl ^ h is the tunnel splitting associated with one path. The explicit forms of 
parameters u and Sj are not relevant for the present calculation; they will be presented in 
subsequent work. 

The upshot of Eq. @ is that the tunnel splitting for flux tunneling can be modulated by 
the application of a gate voltage. When the induced charge is e (half a Cooper pair), the two 
tunneling paths interfere completely destructively and flux tunneling is wholly suppressed. 
Furthermore, when e < q mod (4e) < 3e the tunnel splitting becomes negative. This means 
that the ground and excited states interchange roles: if the ground (excited) state can be 
approximated by (\(f>o) + (— ) \<Pi))/V^ when — e < q mod (4e) < e, where |0 O ) an d \(f>i) 
are the distinct fluxoid states connected by the tunneling paths in Fig. |2|, then it becomes 
(|0o) - (+) |0i»A/2 for e < \q\ mod (4e) < 3e. 

In the limit of large "internal" charging energy, Eq 3> Ej, the physics of suppression of 
flux tunneling for q ~ e remains the same as for Eq <C Ej. The quantitative form of the 
tunneling amplitude is, however, quite different. The difference is due to the fact that for 
Ec 3> Ej the system wave function is delocalized in the ^-direction and flux trajectories 
with all values of A# contribute to the tunneling rate. Since the regime of flux tunneling 
requires the "external" charging energy Eq = (2e) 2 /2C to be smaller than Ej, all energies of 
the flux dynamics are then smaller than the energies of the charge dynamics on the central 
island. In this case, when q — e, only the two charge states of the island, n = and n — 1, 
are relevant for the charge dynamics. The Hamiltonian of the system reduces to: 

O 2 f$ - $ ) 2 e(a - e) 

H = ^+ { OT x) - E + cos(7r$/$ K + ^-^a x - E. sin(vr$/$ )a, , (7) 

where E± = (Eji ±Ej 2 )/2 and the Pauli matrices Oi act on the basis of (|0) ± |l))/v2, the 
symmetric and antisymmetric superpositions of the charge states. 

The potential U(Q) for the evolution of the flux $ can be seen to have two branches, 
depending on the charge-space state of the system: For $ < $0/2, the symmetric charge 
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state ( 1 0) + 1 1))/\/2 provides the lower branch of the potential, whereas for <3> > $o/2, the sign 
of cos(7r$/$o) changes and the lower branch of the potential is given by the antisymmetric 
state (|0) - |l))A/2. This means that to avoid suppression of flux tunneling from a state 
localized in the region $ < $o/2 into the state with $ > $o/2, the system should make 
a transition between the symmetric and antisymmetric charge states in the course of flux 
tunneling. From the wave function of these states in the ^-representation: |^ + | oc cos(#/2) 
and j^-l oc sin(6 l /2), one can see that such a transition corresponds to the shift in 9 by ±7r, 
analogous to the similar shift in the limit of large Ej. 

The Hamiltonian (|7]) shows that the coupling between the two branches of the potential 
is provided by the difference E_ of the Josephson coupling energies and deviations of the 
induced charge q from e. In the absence of coupling, transitions between the two potential 
branches and, as a result, flux tunneling, are suppressed. When the coupling is weak (i.e., 
E_ and q — e are small), transitions between the potential branches takes place in the 
vicinity of the degeneracy point $ = $o/2 and can be described in terms of Landau-Zener 
tunneling. The only difference from the standard situation of Landau-Zener tunneling is that 
the transition is not driven by a classical external force but by the flux motion under the 
tunnel barrier and therefore can be described as occurring in imaginary time. The amplitude 
of such an imaginary-time Landau-Zener transition has been found in the context of the 



dynamics of Andreev levels in superconducting quantum point contacts [f24]]. Adapted for 
the present problem, the expression for the transition amplitude is: 



r(A) \\) \ej • 2£ + [£ £]i/2 ' 

Here T denotes the Gamma function and E is the absolute value of the energy of the 
quasistationary flux state at $ < <3> /2 measured relative to the effective top of the potential 
barrier formed by the crossing at $ = $o/2 of the two branches of the potential. The phase 
(f) of the tunneling amplitude coincides with the phase of the coupling between the potential 
branches in the Hamiltonian (|7|) and is given by 

tan</>= — ^— . (9) 

e{q - e) 

Since the overall amplitude of flux tunneling is proportional to w and since w(X) ~ 
(27TA) 1 / 2 for A < 1, Eq. (^) shows that, when E_ is small, flux tunneling is suppressed 
at q ~ e. Equation (Q) also shows that, similarly to the limit of large Ej, the sign of the 
tunneling amplitude at E- — > changes abruptly when q moves through the point q = e. 
The absolute value of amplitude w (|J) is shown in Fig. |3| as a function of charge q for 
several values of E_. The curves in Fig. |3] indicate that the suppression of the flux tunneling 
amplitude at q ~ e remains pronounced for quite large degree of asymmetry of the Josephson 
coupling energies of the two junctions. 

The suppression of tunneling occurs for arbitrary Ej/Ec, as can be seen from the fol- 
lowing argument. The ground-state degeneracy responsible for the suppression of tun- 
neling is due to two symmetries of Hamiltonian ([I]) that occur for Eji = Ej 2 when 
<& x = $o/2 and q = e. The Josephson coupling terms in the Hamiltonian can be 
written as -2Ej cos(7r$/$ ) cos# = —Ej cos(7r$/$ ) E(l^) (n + 1| + |n + l)(n|). Then 
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the Hamiltonian is symmetric with respect to the following two transformations: (I) 
|n) — > (-l) n \n) , $ — > $ - which requires that $2 = $ /2, and (II) n — > 1 - n, 
which requires that q = e. We now show that one cannot construct an energy eigenfunction 
that is simultaneously an eigenfunction of both of these transformations and, therefore, that 
each energy eigenfunction must be degenerate with another one with which it is related by 
the application of both transformations. Let us write the energy eigenfunction in $, n space: 
|\& (n, $)) = E^n ($) |n). Transformation (I) gives ■?/>„ ($) = ( -l) n ^ n ($ - $), i.e. ^ n ($) 

is an odd (even) function of $ — $o/2 if w is odd (even). Transformation (II) requires that 
ip n ($) = ^i-n (<&), but this is impossible because it equates an even function with an odd 
one. Thus, at <& x = $o/2 and q = e, the eigenfunctions of Hamiltonian (|l|) are all two-fold 
degenerate for arbitrary Ej/Ec ratio. 

Our results are a variant of the Aharonov-Casher effect applied to superconductors, like 
that studied by Reznik and Aharonov ||. It is important to note that in our calculation, 
the relevant charge is not n, which may not even be a good quantum number, but q, the 
charge induced by the gate voltage. Another interesting feature of our results is that the 
amount of flux that tunnels need not be a multiple of $o> as noted above and illustrated in 
Fig. Thus, unlike the case studied in || where ffuxons of flux $ move around a localized, 
external charge, the modulation of flux tunneling in the SQUID does not require that the 
flux $ tunnel between quantized values or that the charge on the island be localized. 

We wish to thank J. Mannik, E. Chudnovsky and A. Guillaume for useful discus- 
sions. This work was supported in part by the NSA and ARDA under ARO contract 
DAAD199910341. 
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FIG. 1. Schematic of the proposed device, an rf-SQUID with the single Josephson junction 
replaced by a Bloch transistor. The charge on the superconducting island of the transistor can be 
induced by a gate voltage. 




FIG. 2. Two-dimensional potential for the system described by Eq. (|T|) with Ej\ = Ej2 = Ej 
and Ej = 0.507 < I ) q/2L. When q = e, destructive interference between the two paths shown leads 
to suppression of flux tunneling. 
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FIG. 3. Absolute value of the transition amplitude between the two branches of the flux poten- 
tial, which is proportional to the amplitude of flux tunneling, as a function of the induced charge 
q at q ~ e. Note that the scale qo = Cy,(2E + ) 1 ^ 2 (EqE) 1 ^ 4: /e of variations of q is much smaller 
than e. From bottom to top, the curves correspond to increasing difference between the Josephson 
energies of the two junctions: E-/{2E + ) 1 l 2 {E Q E) l / i = 0.0 , 0.1 , 0.3. 
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